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Abstract. We report on the spectral properties of a sample of 90 hard X-ray selected serendipitous sources detected 
in 12 XMM-Newton observations with 1 ^ -F2-10 ^ 80 x 10~^* erg cm~'^ s~^. Approximately 40% of the sources 
are optically identified with 0.1 ^ z ^ 2 and most of them are classified as broad line AGNs. A simple model 
consisting of power law modified by Galactic absorption offers an acceptable fit to ~ 65% of the source spectra. 
This fit yields an average photon index of (F) ~ 1.55 over the whole sample. We also find that the mean slope 
of the QSOs in our sample turns out to remain nearly constant ((F) ~ 1.8-1.9) between ^ z ^ 2, with no 
hints of particular trends emerging along z. An additional cold absorption component with 10^^ Nh ^ W^^ 
cm~'^ is required in ~ 30% of the sources. Considering only subsamples that are complete in fiux, we find that 



the observed fraction of absorbed sources (i.e. with Nh 



30%, with little evolution in the range 



2 ^ -F2-10 ^ 80 X 10~ erg cm~ s~ . Interestingly, this value is a factor ~ 2 lower than predicted by the 
synthesis models of the CXB. This finding, detected for the first time in this survey, therefore suggests that most 
of the heavily obscured objects which make up the bulk of the CXB will be found at lower fluxes (_F2-io < 10~^^ 
erg cm~^ s~^). This mismatch together with other recent observational evidences which contrast with CXB model 
predictions suggest that one (or more) of the assumptions usually included in these models need to be revised. 



Key words. Galaxies: active - quasars: general - X-rays: 
galaxies - X-rays: diffuse radiation 

1. Introduction 

Deep pencil-beam observations performed with the new 
generation X-ray telescopes, Chandra (Mushotzky et al. 
2000; Brandt et al. 2001; Rosati et al. 2002) and XMM- 
Newton (Hasinger et al. 2001), have resolved the bulk of 
the hard (2-10 keV) cosmic X-ray background (CXB) into 
the integrated contribution from discrete sources, pushing 
the detection limit down to values which are several orders 
of magnitudes fainter than in previous surveys (Moretti 
et al. 2003). Moreover, their excellent angular resolutions 
have allowed the unambiguos identifications of most X- 
ray sources with active galactic nuclei (AGNs) (Hasinger 
2003) providing a unique tool to investigate in detail the 
formation and evolution of accretion powered sources (i.e. 
black holes) over cosmic time as well as the physical con- 
nections between nuclear activity and the host galaxy (Silk 
& Rees 1998; Franceschini et al. 1999). 

Stacked spectra of these faint sources exhibit a flat 
slope in full agreement with the spectral shape of the un- 
resolved CXB in the 2-10 keV band (F = 1.4; Gendreau 
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et al. 1995), thus resolving the so-called "spectral para- 
dox" (De Zotti et al. 1982). Surprisingly, the observed red- 
shift distribution of these sources seems to peak at z < 1 
(GiUi 2003, Hasinger 2003), at odds with the expectations 
of most popular CXB synthesis models (Comastri et al. 
1995; Gilli et al. 2001) which predict a peak at z ^ 1.5-2. 

Since the diffuse CXB emission is now definitively re- 
solved into point sources, the general interest has moved to 
accurately constrain the physical and evolutionary prop- 
erties of the different classes of X-ray sources. However, 
these extremely deep pencil-beam exposures detect the 
majority of the sources with a very poor counting statistics 
which prevents from inferring accurate object-by-object 
X-ray spectral properties: such information is nevertheless 
essential to reveal the physical conditions and the geome- 
try of the matter in the circumnuclear region and, hence, 
how AGNs ultimately work. 

We have therefore undertaken a wide-area search 
(Piconcelli et al. 2002; hereafter Paper I) aimed at ex- 
ploring the individual spectral properties of moderate to 
faint hard X-ray selected sources. The best way to achieve 
this goal is to exploit the large collecting area of XMM- 
Newton imaging detectors in order to collect the largest 
number of targets with X-ray spectra of good quality. 
It is worth noting that our survey samples a 2-10 keV 
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Table 1. Journal of the XMM-Newton observations. 



No. 


Field 


R.A. 


Dec. 


Date 


Obs. ID 


Exposure ( 




Filter^ 
















Mi 


M2 


I\\ Mi M2 


1 


PKS0312-770'' 


03 11 55.0 


-76 51 52 


2000-03-31 


0122520201 


26000 


25000 


24000 


To To Tc 


2 


MS1229.2+6430" 


12 31 32.0 


+64 14 21 


2000-05-20 


0124900101 


22900 


18600 


22900 


Th Th Th 


3 


IRAS13349+2438'' 


13 37 19.0 


+24 23 03 


2000-06-20 


0096010101 


— 


41300 


38600 


- MeTh 


4 


Abel! 2690° 


00 00 30.0 


-25 07 30 


2000-06-01 


0125310101 


21000 


16600 


15300 


Me Me Me 


5 


MS 0737.9+744" 


07 44 04.5 


+74 33 49 


2000-04-12 


0123100101 


15000 


17800 


26100 


Th Th Th 


6 


Markarian 205" 


12 21 44.0 


+75 18 37 


2000-05-07 


0124110101 


17000 




14800 


Me - Me 


7 


Abell 1835" 


14 01 02.0 


+02 52 41 


2000-06-27 


0098010101 


22900 


23700 


26400 


Th Th Th 


8 


PHL5200'' 


22 28 30.6 


-05 18 32 


2001-05-28 


0100440101 


40200 


34900 


35000 


Tc Tc Tc 


9 


SDSSJ1044-0125'' 


10 44 31.8 


-01 25 09 


2000-05-21 


0125300101 


37500 


28900 


28800 


Th Th Th 


10 


Lockman Hole'' 


10 52 43.0 


+57 28 48 


2000-04-27 


0123700101 


33500 


27000 


29900 


Th Tc Th 


11 


NGC253'' 


00 47 37.1 


-25 17 41 


2000-06-03 


0125960101 


34100 


31600 


30300 


Me Me Th 


12 


LBQS2212-1759'' 


22 15 30.9 


-17 44 14 


2001-11-17 


0106660601 


80500 






Th - - 



'Optical blocking filters used during observations: Th=thin, Me=medium and Tc=thick. 
References: " Paper I; this work. 



flux range {F2-10 ^ 10~^^+10^^'' erg cm^^ s^^) which 
is mostly uncovered by the ultradeep surveys designed to 
resolve the CXB at much fainter flux levels. In Paper I 
we presented the first results of an initial sample of 41 
serendipitous X-ray sources selected from seven XMM- 
Newton observations with moderate (~ 20-40 ks) expo- 
sures. 

In this work, thanks to the addition of further five 
deeper XMM-Newton observations (up to ~ 80 ks), we 
extend our detection limit down to F2_io '--^ 1 x 10^^^ 
erg cm~^ s^^. In this way, we have been reached values 
closer to the knee of the hard logN-logS distribution, i.e. 
where most of the sources accounting for the bulk of the 
CXB reside (Moretti et al. 2003). Accordingly, we are also 
able to place stronger constraints on some input /output 
parameters of synthesis models of the CXB than in Paper 

1. In particular, we are able to provide here a sounder 
estimate of the fraction of absorbed objects in the faint 
hard X-ray population. 

2. XMM-Newton observations and data reduction 

The present study is based on a set of 12 observations 
carried out by the XMM-Newton satellite (Jansen et 
al. 2001) from March 2000 to November 2001, including 
Performance/ Verification phase. Target of Opportunity 
and Guest Program observations available in the XMM 
Newton Science Archive. The imaging and spectroscopic 
measurements are taken by the European Photon Imaging 
Camera {EPIC) consisting of one PN back-illuminated 
CCD array (Struder at al. 2001) and two MOS front- 
illuminated CCD arrays (Turner et al. 2001). These fields 
are chosen for their high galactic latitudes ( |&| ^ 30 
deg) and exposures ( ^ 15 ks) which enable us to collect 
a large sample of cosmic hard X-ray serendipitous sources 
without heavy contamination from our Galaxy. 

Table 1 lists the name and the coordinates of the tar- 
get sources, together with the epoch and the identifica- 



tion number for each of the 12 XMM Newton observation. 
Results from the analysis of the first seven observations 
listed in Table 1 (i.e. from No. 1 to No. 7) are reported in 
Paper I. 

For the five newly included EPIC observations (i.e. 
observations from No. 8 to No. 12 in Table 1) we have ap- 
plied the same data reduction and detection procedures 
described in Paper I. Similarly, we have included in the 
final sample only those 49 hard X-ray selected serendip- 
itous sources which satisfy the same selection criterion 
as used in Paper I (which provides at least ^ 100 net 
counts in the 2-10 keV band once all the three EPIC cam- 
eras are taken into account). The final complete catalog 
of the 90 X-ray sources is listed in Table 2 together with 
their XMM-Newton coordinates. Approximately 40% of 
the sources (i.e. 37 out of 90) are optically identified with 

0. 1 ^ z ^ 2 and most of them are classified as broad line 
AGNs. 

We further create from this basic sample two subsam- 
ples: the BRIGHT sample which includes 42 X-ray sources 
with! F2-W > 5 X 10"" erg cm'^ s^^ and the FAINT 
sample with 22 sources having 2 < F2-10 < 5 x lO"^** erg 
cm~^ s~^. This latter sample has been obtained taking 
into account only the last 5 fields of Table 1, i.e. those 
with the longest exposures, for which we could estimate a 
flux hmit of F2-10 ~ 2 X 10~^^ erg cm~^ s~^. This flux 
limit has been calculated from the shortest observation, 

1. e. the Lockman Hole with an exposure of 33.5 ks. By do- 
ing so, also the FAINT sample is therefore complete down 
to this flux limit. 



^ Fluxes of the sources are calculated using the best fit spec- 
tral model (see Paper I and Table 6). 
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Table 2. The hard X-ray selected sample. 



N 


Source name 


R.A. 


Declination 


z 


Classification 


7?(t) 








(J2000) 


(J2000) 






(mag) 


(mjy) 






PKS 0312-770 field 












1 


CXOUJ031015.9- 765131 


03 10 15.3 


-76 51 




1.187 


BL AGN" 


17.6 




2 


CXOUJ031209. 2- 765213 


03 12 08.7 


-76 52 




0.89 


BL AGN° 


18.2 


_ 


3 


CXOUJ031238.9-765134 


03 12 38.8 


-76 51 


31 


0.159 


Galaxy" 


17.7 


_ 


4 


CXOUJ031253.8-765415 


03 12 53.5 


-76 54 


1 


0.683 


Red QSO" 


22.0 


_ 


5 


CXOUJ031312.1-765431 


03 13 11.5 


-76 54 




1.124 


BL AGN" 


18.3 


_ 


6 


CXOUJ031314.5- 765557 


03 13 14.2 


-76 55 


54 


0.42 


BL AGN" 


19.1 


_ 


7 


XMMUJ03091 1 .9-765824 


03 09 11.6 


-76 58 


24 


0.268 


Sey 2^ 


19.1 




8 


XMMUJ031049.6- 763901 


03 10 49.5 


-76 39 


m 


0.380 


BL AGN*" 


18.6 




9 


XMMUJ031105. 1-765156 


03 11 05.1 


-76 51 


oo 




No cl. 










MS1229.2+6430 field 












10 


XMMUJ123110.6-I-641851 


12 31 10.6 


+64 18 


Ol 




No cl. 


18.7 




11 


XMMUJ123116. 3-^641114 


12 31 16.3 


+64 11 


1 A 




No cl. 






12 


XMMUJ123218. 6-^640309 


12 32 18.6 


+64 03 


HQ 




No cl. 


20.0 




13 


XMMUJ123214.2-F640459 


12 32 14.2 


+64 04 






No cl. 






14 


XMMUJ123013.4+642505 


12 30 13.4 


+64 25 


UD 




No cl. 


15.9 




15 


XMMUJ123049. 9-1-640845 


12 30 49.9 


+64 08 






No cl. 


18.6 




16 


XMMUJ123058. 54^641726 


12 30 58.5 


+64 17 






No cl. 


20.0 






IRAS13349-F2438 field 












17 


XMMUJ133730.8-F242305 


13 37 30.8 


+24 23 


05 


_ 


No cl. 


19.39 


_ 


18 


XMMUJ133649.3+242004 


13 36 49.3 


+24 20 


04 





No cl. 


19.97 





19 


XMMUJ133807.4-F242411 


13 38 07.4 


+24 24 


11 





No cl. 


18.17 





20 


XMMUJ133747.4-F242728 


13 37 47.4 


+24 27 


28 


— 


No cl. 


19.5 


— 


21 


XMMUJ133712.6+243252 


13 37 12.6 


+24 32 


52 


— 


No cl. 




— 






Abell 2690 


field 












22 


XMMUJ000031.7-255459 


00 00 31.7 


-24 54 


59 


0.283 


BL AGN*- 


17.7 


— 


23 


XMMUJ000122.8-250019 


00 01 22.8 


-25 00 


19 


0.968 


BL AGN*" 


18.7 


69.2^ 


24 


XMMUJ000027.7-250441 


00 00 27.7 


-25 04 


41 


0.335 


BL AGN*" 


18.6 





25 


XMMUJOOOlOO.0-250459 


00 01 00.0 


-25 04 


oy 


0.851 


BL AGN*" 


21.9 


130^*' 


26 


XMMUJ000102. 5-245847 


00 01 02.5 


-24 58 




0.433 


BL AGN*" 


20.3 




27 


XMMUJ000106.8-250845 


00 01 06.8 


-25 08 
















MS 0737.9+744 field 












28 


1E0737.0-F7436 


07 43 12.5 


+74 29 


35 


0.332 


BL AGN= 


16.4 




29 


XMMUJ074350.5+743839 


07 43 50.5 


+74 38 


39 




No cl. 







30 


ISAX J0741.9-h7427 


07 42 02.2 


+74 26 


24 





No cl. 


19.0 





31 


XMMUJ074351.5-F 744257 


07 43 51.5 


+74 42 


57 





No cl 


20.0 





32 


XMMUJ074401.5+743041 


07 44 01.5 


+74 30 


41 


_ 


No cl. 


_ 


87 






Markarian 205 field 












33 


MS1219.9+7542 


12 22 06.6 


+75 26 


14 


0.238 


NELG'' 


16.91 


_ 


34 


MS1218.6-f7522 


12 20 52.0 


+75 05 


29 


0.646 


BL AGN'* 


17.7 


_ 


35 


XMMUJ122258.3-F 751934 


12 22 58.3 


+75 19 


34 


0.257 


NELG"* 




_ 


36 


XMMUJ122351.3+752224 


12 23 51.3 


+75 22 


24 


0.565 


BL AGN'^ 




_ 


37 


NGC4291 


12 20 15.9 


+75 22 


09 


0.0058 


Galaxy'* 


1 1 .7 









Abell 1835 field 












38 


XMMUJ140127.7-h025603 


14 01 27.7 


+02 56 


03 


0.265 


BL AGN*" 


19.7 


I.54W 


39 


XMMUJ140053.0+030103 


14 00 53.0 


+03 01 


uo 


0.573 


BL AGN*" 


19.7 




40 


XMMU J 140130. 7+024529 


14 01 30.7 


+02 45 


9Q 




No cl. 






41 


XMMUJ140145.0+025330 


14 01 45.0 


+02 53 


oU 


_t 


Galaxy''' + 


17.9 








PHL 5200 


field 












42 


XMMUJ222814.0-051621 


22 28 14.0 


-05 16 


21 




No cl. 


19.73 




43 


XMMUJ222814.9-052418 


22 28 14.9 


-05 24 


18 





No cl. 







44 


XMMUJ222834.5-052150 


22 28 34.5 


-05 21 


50 




No cl. 






45 


XMMUJ222822. 1-052732 


22 28 22.1 


-05 27 


32 




No cl. 






46 


XMMUJ222850.5-051658 


22 28 50.5 


-05 16 


58 




No cl. 






47 


XMMUJ222823.6-051308 


22 28 23.6 


-05 13 


08 




No cl. 






48 


XMMUJ222905. 2-051432 


22 29 05.2 


-05 14 


32 




No cl. 






49 


XMMUJ222732. 2-051644 


22 27 32.2 


-05 16 


44 




No cl. 


19.92 


3.1 



4 
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Table 2. continued 



N 


Source name 


R.A. 


Declination 


z 


Classification 










(.J2()()()) 


(.1201)1)) 






mag 


m.Jy 


50 


XMMUJ222826.7-051821 


22 28 26.7 


-05 18 21 


— 


No el. 


— 


134 




SDSS J1044-0125 field 










51 


XMMUJ104451.3-012229 


10 44 51.3 


-01 22 29 


- 


No cl. 


- 


- 


52 


XMMUJ104445. 1-012420 


10 44 45.1 


-01 24 20 


— 


No cl. 


— 


— 


53 


2QZJ104424.8-013520 


10 44 25.1 


-01 35 20 


1.57 


BL AGN 


18.61 


— 


54 


XMMUJ104509.4-012441 


10 45 09.4 


-01 24 41 


- 


No cl. 


19.6 


- 


55 


XMMUJ104456.0-012533 


10 44 56.0 


-01 25 33 


- 


No cl. 


19.0 


- 


56 


XMMUJ104441.9-012655 


10 44 41.9 


-01 26 55 


- 


No cl. 


- 


- 


57 


2QZJ104522. 0-012845 


10 45 22.3 


-01 28 54 


0.782 


BL AGN 


18.5 


- 


58 


XMMUJ104444.6-013315 


10 44 44.6 


-01 33 15 


— 


No cl. 


18.4 


3.74 




The 


Lockman Hole field 










59 


RXJ105421.1-F572545 


10 54 21.1 


+57 25 45 


0.205 


Sey 1.9 (-^^ 


18.3 


0.8 


60 


RXJ105316.8-h573552 


10 53 16.8 


+57 35 52 


1.204 


BL AGN 


19.0 


0.26* 


61 


RXJ105239.7-^572432 


10 52 39.7 


+57 24 32 


1.113 


BL AGN 


18.0 


0.14 


62 


RXJ105335.1+572542 


10 53 35.1 


+57 25 42 


0.784 


BL AGN 


19.78 


— 


63 


7C 1048+5749 


10 51 48.8 


+57 32 48 


0.99 


NL AGN 


22.9 


15.39* 


64 


RXJ105339.7+573105 


10 53 39.7 


+57 31 05 


0.586 


BL AGN '•f 


19.4 


- 


65 


XMMUJ105237.8-I-573322 


10 52 37.8 


+57 33 22 


0.707 


NL AGN ^f'> 


22.6 


59.45* 


66 


RXJ105331.8+572454 


10 53 31.8 


+57 24 54 


1.956 


BL AGN 


19.99 


- 


67 


RXJ105350.3+572709 


10 53 50.3 


+57 27 09 


1.720 


BL AGN 


20.15 


- 






NGC 253 field 










68 


RXJ004759.9- 250951 


00 47 59.9 


-25 09 51 


0.664 


BL AGN 


17.4 


- 


69 


XMMUJ004722. 5-251202 


00 47 22.5 


-25 12 02 


- 


No cl. 


- 


— 


70 


RXJ004722. 9-251053 


00 47 22.9 


-25 10 53 


1.25 


BL AGN ' 


18.1 


— 


71 


RXJ004647.2-252152 


00 46 47.2 


-25 21 52 


1.022 


BL AGN '•<>'> 


20.17 


_ 


72 


XMMUJ004818.9-251505 


00 48 18.9 


-25 15 05 


_ 


No cl. 


_ 


_ 




LBQS 2212-1759 field 










73 


XMMUJ221536. 5-173357 


22 15 36.5 


-17 33 57 


_ 


No cl. 


18.0 


_ 


74 


XMMUJ221510. 7-173644 


22 15 10.7 


-17 36 44 


_ 


No cl. 


_ 


_ 


75 


XMMUJ221604.9-175217 


22 16 04.9 


-17 52 17 


_ 


No cl. 


20.42 


_ 


76 


XMMUJ221557.8-174854 


22 15 57.8 


-17 48 54 


_ 


No cl. 


_ 


_ 


77 


XMMUJ221623. 1-174055 


22 16 23.1 


-17 40 55 


_ 


No cl. 


_ 


13.6 


78 


XMMUJ221519.4-175123 


22 15 19.4 


-17 51 23 


— 


No cl. 


— 


— 


79 


XMMUJ221453.0-174233 


22 14 53.0 


-17 42 33 


— 


No cl. 


— 


— 


80 


XMMUJ221518.8-174005 


22 15 18.8 


-17 40 05 


— 


No cl. 


— 


— 


81 


XMMUJ221602.9-174314 


22 16 02.9 


-17 43 14 


— 


No cl. 


— 


— 


82 


XMMUJ221623. 7-174722 


22 16 23.7 


-17 47 25 


— 


No cl. 


20.99 


— 


83 


XMMUJ221602. 9-174314 


22 16 02.9 


-17 43 14 


— 


No cl. 


— 


— 


84 


XMMUJ221537.6-173804 


22 15 37.6 


-17 38 04 


- 


No cl. 


- 


- 


85 


LBQS 2212-1747 


22 15 15.0 


-17 32 24 


1.159 


BL AGN 


17.3 




86 


XMMUJ221623. 5-174317 


22 16 23.5 


-17 43 17 




No cl. 


20.9 




87 


XMMUJ221550.4-175209 


22 15 50.4 


-17 52 09 




No cl. 


18.1 




88 


XMMUJ221533.0-174533 


22 15 33.0 


-17 45 33 




No cl. 






89 


XMMUJ221456.7-175054 


22 14 56.7 


-17 50 54 




No cl. 






90 


XMMUJ221523.7-174323 


22 15 23.7 


-17 43 23 




No cl. 


20.44 





Optical classifications and redshifts are taken from: '■'^^ Fiore et al. (2000), Fiore et al. 2003 (F03; in preparation), Wei 
et al. (1999), (''^ AXIS (e.g. Barcons et al. 2001), Veron-Cetty & Veron (2001), '■f^ Mainieri et al. (2002), Vogler & 
Pietsch (1999). (f) There are two possible candidates for the identification of this sources: an elliptical galaxy at z = 0.251 or 
an elliptical galaxy at z =0.254 (F03). Magnitude in the R band. Photometric data are taken from the USNO catalog or 
F03 whenever available. ^-^^ Flux density at 1.4 GHz (i.e. 20 cm). Data are taken from FIRST and NVSS on-line catalogs. '■*^ 
Radio loud (RL) object. 



3. Spectral analysis 

In this Section we focus on the spectral analysis of the 
X-ray sources selected in the five new XMM-Newton ob- 



servations i.e. sources from No. 42 to No. 90 in Table 2. 

Detailed results about the analysis of the first 41 sources 
listed in Table 2 can be found in Paper L We have per- 
formed the spectral analysis in the 0.3-10 keV band, 
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choosing the background region in the same detector chip 
and with the same extraction radius of the source region. 
Wc adopt as standard a circular extraction region of 35 
arcsecs radius, both for PN and MOS, shortened if other 
X-ray sources or CCD gaps are present inside this region. 
If during the observation the optical filter of the two MOS 
cameras are the same we combine together their spectra. 
Whenever both datasets are available, a joint spectral fit- 
ting using the data of both PN and MOS is carried out. 
The XSPEC v. 11.0.1 software package has been used to 
analyse all the background-subtracted source spectra. In 
order to permit fitting, we use a minimum spectral 
group size of 20 events per data points. However, in the 
case of faint sources with < 400 counts in the broad band 
0.3-10 keV, wc rcbin the data so that there are at least 15 
counts in each bin and we applied the Gehrels weighting 
function in the calculation of (Gehrels 1986) since it 
is a better approximation in the calculation of when 
the number of net counts is small. For the spectral analy- 
sis we have used the latest known response matrices and 
calibration files (January 2002) released by the XMM~ 
Newton Science Operations Centre, taking into account 
the type of optical filters applied at the top of the tele- 
scopes during the observations (sec Table 1). 

Throughout this paper we adopt Hq = 50 km s~^ 
Mpc~^ and Qo = Q for the calculation of the luminosi- 
ties. Unless stated otherwise, the errors refer to the 90% 
confidence level for one interesting parameter (i.e. Ax^ = 
2.71; Avni 1976). 

3.1. Spectral fitting 
3.1.1. Basic models 

We begin the spectral analysis by fitting the spectra with 
a simple power law plus Galactic absorption (SPL) model. 
This basic spectral parameterization allows us to look for 
any evidence of absorption and/or excess emission fea- 
tures: it also provides useful indications about the mean 
slope of the continuum at these hard X-ray flux levels 
(see Sect. 6). Results of these fits are displayed in Tabic 3, 
while the spectrum of each source of the sample together 
with the relative data-to-model ratios can be found in 
PiconcelH (2003). 

We find that the values of x^^^ are statistically accept- 
able in most cases, thus suggesting that the SPL model 
provides a reasonable description of the data for the ma- 
jority (29 out of 49, i.e. 60%) of the sources. It must be 
borne in mind that for the faintest objects, the low qual- 
ity data prevent a very detailed modeling of some features 
(i.e. lines, edges) which are possibly present in their spec- 
tra: as consequence the SPL model provides a good de- 
scription of the overall spectral shape even if it is not the 
most appropriate. 

Some sources are, however, clearly not satisfactorily 
fitted by the SPL model because either the data-to-model 
ratio residuals are present or y^^^ » 1 (see Table 3). 
Furthermore, spectra with flat photon index (F ^ 1.3- 



Table 3. Spectral fitting results - I. Fits with a single power 
law plus Galactic absorption model (SPL). 



N 


r 


X;'ed/(d.0.f.) 




Best-fit 






PHL 5200 field (Ng"-'= 5.3 


X 


10^° cm-^) 




42 


1 -79+0.16 
' ^-0.16 


63/f42"l 




Yes 




43 


, o 7+0. 44 


1 03/(26) 




No 




44 




1 54/1^521 




Yes 




45 


l.UU_Q 20 


91/(42) 




No 




46 


1 7r+0.22 
-■■•'^-0.22 


66/(32) 




Yes 




47 

48 


i.84«:- 
i.oot°:l° 


96/1^651 

1 45/(74) 

±.'±0/ y 1 'i) 




Yes 

No 




49 


1-841H^ 


1.05/(191) 




Yes 




50 


0.85lH^7 


1.35/(176) 




No 




SDSS J1044-0125 


field (Ng=''= 


:4.2 X 10^" cm 


"') 


51 


1 90+0.10 
J--^**-o.io 


1.01/(85) 




No 




52 


1 c 0+0. 13 
J^-^'J-0.13 


0.99/(44) 




Yes 




53 


-■■•"^-0.10 


0.93/(86) 




Yes 




54 


2.05lg:- 


0.84/(88) 




Yes 




55 


2.03l°:« 


0.85/(37) 




Yes 




56 


U.4Z_Q 24 


2.10/(19) 




No 




57 


2 04+°-^^ 


0.73/(100) 




Yes 




58 


O.DD_o.28 


1.46/(28) 




No 




Lockman Hole field (Ng'''= 5.5 


X 10^^ cm- 


') 


59 


1 91 +0.08 
^•^^-0.08 


1.54/(361) 




No 




60 


1 70+0.05 
.■■•'^-0.05 


1.16/(178) 




Yes 




61 


2 42+0.10 
^•^^-0.10 


1.28/(126) 




Yes 




62 


9 no+o-os 

^•'-'^-0.08 


1.26/(135) 




Yes 




63 


71 +0-19 


1.10/(38) 




No 




64 


^"J^-O.ll 


1.10/(121) 




Yes 




65 


r, 99+0.40 

-U.ZZ_g 4g 


1.07/(11) 




No 




66 


1 Q7+0.14 

-0.14 


1.16/(61) 




Yes 




67 


1 74+0.28 


0.91/(29) 
[ (Ng"'= 1.5 




Yes 






NGC 253 field 


X 


10^" cm-^^) 




68 


-■■•^"-o.og 


1.11/(119) 




Yes 




69 


-0 24+"-''* 


0.69/(17) 




No 




70 


1 09 + 0. 18 
-l-^^-0.18 


1.00/(69) 




Yes 




71 




1.29/(85) 




Yes 




72 


-o.2ol«:« 


1.07/(14) 




No 




LBQS 2212-1759 field (Ng^'= 


2.4 X 10^° cm 


"') 


73 


2.i4l°:^2^ 


1.10/(63) 




Yes 




74 


-o.i8l°:g 


0.60/(19) 




No 




75 


2.08^°:f2 


1.06/(60) 




Yes 




76 


1 n8+o-44 

-■-1^^-0.36 


0.83/(22) 




No 




77 


71+0-31 

J^-0.31 


0.87/(31) 




No 




78 


2 10+°-" 

^•J^^-0.17 


1.15/(62) 




Yes 




79 


-■■•"^-0.24 


1.10/(32) 




No 




80 


-■■•-1^-0.32 


1.02/(29) 




No 




81 


2 26+"-ii 
'^•^"-o.ii 


1.11/(40) 




Yes 




82 


2 03+"-i^ 


1.08/(133) 




Yes 




83 




0.84/(49) 




Yes 




84 


U.OU„o 28 


1.73/(22) 




No 




85 

86 


9 f;9+0.09 
^•"^-0.10 

i-86^;]:i;: 


0.91/(162) 

0.97/(135) 




Yes 

Yes 




87 

88 


9 99+0.13 
^•^^-0.13 
1 10+0.22 


0.82/(88) 
1.21/(24) 




Yes 
No 




89 


1 06+°" 

l.UO_Q j^4 


1.29/(63) 




No 




90 


9 99+0.11 
^•^^_0.11 


0.88/(58) 




Yes 
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1.4) indicate the likely presence of intrinsic obscuring ma- 
terial which suppresses the soft X-ray continuum. 

We have therefore refitted each spectrum applying a 
power law plus an additional absorption component (in 
source-frame if the rcdshift is known) . This spectral model 
will be indicated as APL hereafter. In Table 4 we report 
the relative spectral parameters for those source spectra 
showing a significant improvement at > 95% confidence 
level according to an i^— test once compared to the SPL 
model fit. Values of the F-statistic and the corresponding 
significance level arc also listed in this Table. 

We also include in Table 4 three faint objects (i.e. No. 
72, No. 74, No. 88) despite the fact that the APL model 
for them is not significantly better than SPL. Indeed their 
extremely flat SPL spectra strongly suggest the presence 
of heavy absorption but owing to the relatively poor statis- 
tics, it is not possible to accurately constrain this compo- 
nent. 

In particular, the application of the APL model to 
source No. 74 reveals strong obscuration (Nh ~ 2 x 10^^ 
cm~^) but the relative spectral index still remains very flat 
and loosely constrained because of a likely simultaneous 
contribution of several unresolved spectral components. 
We therefore fix the photon index F = 1.9, i.e. the mean 
value observed in bright AGNs (Nandra & Pounds 1994), 
to obtain an estimate of the absorption column density 
value in this source (Nh ~ 10^2-23 cm~^, see Table 4). 

Using the APL model we find column densities span- 
ning from ~ 10^^ to ~ 2 X 10^^ cm~^ (see Table 4): in 
particular, broad line objects have low amount of cold 



Table 4. Spectral fitting results - 11. Fits with a single power 
law plus extra absorption component model (APL). 



N 


r 


N„ 

(iO-^ cm -) 


F't'/C.l."' 


Best fit 


45 


1 fji+o.ei 


q c:+4.0 


4.8/96% 


Yes 


48 


1 Q2+°-28 
-l-f^_0.23 


f>A+_\\ 


70/ >99.9% 


Yes 


50 


1 09+°-i° 


1 S+O-9 


15.2/98.5% 


No 


51 


1 04+0-14 


9 90+0.50 
^•^'^-0.50 


49.6/ >99.9% 


Yes 


56 


1 74+0.31 


ii.otii," 


19/ >99.9% 


Yes 


58 


1 10+0.49 
J--io_0.49 




4.5/97.5% 


No 


59 


1 or+0.09 

1.00_Q Q9 


9 1 +0.3 
^■^-0.3 


207/ >99.9% 


Yes 


63 


1 07+0.44 
^-•■J' -0.51 


r,q + 20.7 


5.3/97%. 


No 


65 


1 fi«+0-72 


1 1 S r.+ss-o 

ii8.0_i„g 


4.8/96% 


Yes 


69 




c 19+7.88 


15.7/ >99.9% 


Yes 


72 


1 oq+2.16 

-L.OO 'j'j 


<67.5 


2.0/82% 


No 


74 


1.9f.' 




-/- 


Yes 


76 


1 qq+i-58 


3 3+''-^ 


4.9/97% 


Yes 


79 


q Qq + 1.52 


q 1 +2.0 
•J-1-1.3 


12.4/ >99.9% 


No 


84 




10.21^^2* 


7.6/98.5% 


Yes 


88 


1 c:n+0.46 
X.OU_Q 4g 




2.70/90% 


Yes 


89 


1 Q1+0.48 
^•^^-0.28 


q 9+1.5 


19.4/ >99% 


Yes 



' F-statistic value. * Confidence level with respect to model 
SPL (see Table 3) using the F-statistic. 



absorption (i.e., Nr < 10^^ cm ^), similarly to what 
found in Paper I. Interestingly both 7C 1048+5749 and 
XMMU J105237.8+573322 (source No. 63 and No. 65, re- 
spectively) show X-ray luminosities exceeding 10^** erg/s 
as well as a column density Nh > 10^^ cm~^, thus becom- 
ing candidates to be type 2 QSOs (e.g. Sect. 5, Mainieri 
et al. 2002). Note that the values of spectral parameters 
r and Nh derived by our analysis of the sources in the 
Lockman Hole field (i.e. sources from No. 59 to No. 67) 
fully agree with those obtained by Mainieri et al. (2002) 
using a longer (~ 100 ks) XMM-Newton observation. 

As expected the introduction of an additional absorp- 
tion component produces a significant steepening of the 
continuum slope in most of the sources listed in Table 4. 
However, a sizeable number of objects (5 out of 17) still 
have fiat spectra with F ^ 1.2 and/or exhibit residuals 
in their data-to-model ratios. Thus a further and more 
detailed analysis has been carried out in order to take 
into account also these additional spectral features (see 
Sect. 3.1.2). 

Finally, as already done in Paper I, we also fit all the 
spectra with the APL model fixing F = 1.9 (and ^; = 1 for 
the optically unidentified sources) in order to overcome a 
possible underestimation of the intrinsic column densities 
in those sources with the lowest statistics and/or with- 
out redshift information^. We choose z = 1 on the basis 
of the findings reported in recent optical follow-ups of ul- 
tradeep X-ray surveys (Hasinger et al. 2003; Cowie et al. 
2003) which suggest a peak at z ;S 1 in the redshift dis- 
tribution of the sources making the CXB. Results of this 
spectral fitting will be discussed in Sect. 6 in the frame 
of the observational constraints on the predictions of the 
synthesis model of the CXB. 

3.1.2. More complex models and peculiar sources 

Although intrinsic absorption suppresses a sizeable frac- 
tion of the soft X ray primary continuum, many Type 2 
AGNs are characterized by a soft-excess component which 
is either originating in a circumnuclear diffuse starburst 
and/or is due to reprocessed emission scattered along our 
line of sight by a photoionized gas located just above the 
obscuring torus (Turner et al. 1997, Kinkhabwala et al. 
2002) . This is the case for seven X-ray sources in our sam- 
ple (i.e. Nos. 50, 58, 63, 72, 77^, 79'' and 80; see Table 4), 

^ The effective column density N^^-'^ has the following red- 
shift dependance: N^^-^ oc Nh(1 + zf-'^ (Barger et al. 2002). 

Source No. 77 has been included hero due to its very flat 
photon index derived by the SPL model. Accordingly, this 
source is likely obscured by a large amount of absorption 
but the data quality does not allow a more accurate spectral 
modeling. 

* Although the statistical improvement is not so significant 
source No. 79 has been included here because after the intro- 
duction of a soft-excess component the resulting photon index 
is r « 2 (i.e. a value commonly found in AGNs) instead of the 
unusual steep slope derived by the APL model (F w 3) 
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for which we have therefore included such a component 
(significant at ^ 90% confidence level) in their best fit 
model. A thermal Raymond-Smith component (labelled 
with TM in Table 5) is required in source Nos. 79 and 
80; while an additional power law (labelled with PL in 
Table 5) is added for the latter five remaining objects. The 
metallicity of the thermal component is fixed to the solar 
value while the spectral index of the second power law is 
put equal to the value found for the hard X-ray primary 
power law, as expected in the case of a scattered com- 
ponent (Turner et al. 1997). An example of an absorbed 
source (No. 50) for which we have applied an additional 
power law spectral parameterization is shown in Figure 1. 

The resulting average increase of the intrinsic column 

density value due to the addition of a soft excess compo- 
nent is (Nh) 




channel energy (keV) 

Fig. 1. The EPIC {PN and MOS) spectrum of the unidenti- 
fied source No. 50 in the PHL 5200 FOV with the soft-excess 
component fitted with a scattered power law model (PL in 
Table 5). 



Finally, a clear warm absorber signature is present in 
the SPL spectrum of the unidentified source No. 43: we 
have therefore added in its fitting model an absorption 
edge to parameterize this feature. The improvement in the 

is significant at ^ 96% confidence level with a resulting 
observed-frame energy^ for the edge Egdge = 0.55tQ;o7 
keV, likely due to highly ionized OVII/OVIII as commonly 
found in many Seyfert Is (Reynolds 1997). 

In Table 6 flux in the 0.5-2 keV and 2-10 keV band, 
2-10 keV luminosity and best fit model are listed for all 
sources presented in this Section. 



4. Results on the whole sample 

Adding the five XMM-Newton exposures presented in 
Sect. 3 (i.e. observations Nos. 8 to 12 in Table 1), the 
number of hard X-ray selected sources in our sample has 
increased from 41 (in Paper I) to 90. In this Section we 
present the results obtained by taking into account the 
whole sample as listed in Table 2. The present XMM- 
Newton observations yield the first 0.3-10 keV spectrum 
of a large fraction of the X-ray sources investigated in 
this work, since most of them have not been detected by 
previous less sensitive X-ray telescopes. The results in- 
ferred by the spectral analysis of the entire sample can 
be briefiy summarized as follows. For about 65% (i.e. 60 
out of 90) of the X-ray sources the SPL model represents 
an acceptable description of their spectra. 26 (out of 90) 
sources require the introduction of a significant absorption 
component (see Table 4, Table 5 and Paper I). The result- 
ing Nh distribution is shown in Figure 2. Furthermore, 13 
sources require a more complex fitting model than SPL or 
APL to account for a soft excess component (11 out of 13) 
or for the presence of warm absorber signatures (2 out of 
13). 

Measured values of the hard X-ray flux range from 

with more than 50 

< 5x IQ-^^ erg 
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~ 1 to ~ 80 X 10"^'' erg cm~^ s~ 
out of 90 sources (i.e. 55%) at F2- 
cm~^ i.e. flux levels almost unexplored by the X-ray 
telescopes operating before XMM-Newton. As expected 
on the basis of our selection criterion, in the soft X-ray 
band we detect sources in a broader flux range i.e. from 
^0.5-2 « 70 down to « 0.04 x 10"" erg cm-^ 8"^ The 
absorption-corrected 2-10 keV luminosities span from w 
2 X 10^° erg s~^ to !=s 5 x 10'*^ erg s~^ in agreement with 
the optical classification of the identified sources in the 
sample. All but one (i.e. NGC 4291, n. 37) sources have a 
^2-10 ^ iC^ erg s~^ typical of AGN: the two optically 



Log Column density cm — 2 



^ Assuming that the observed edge is due to OVII(OVIII), 
we infer a redshift z = 0.35(0.81) for this X-ray source. 



Fig. 2. Distribution of the absorption column densities for the 
26 sources in the sample requiring the introduction of a signif- 
icant absorption component (see Table 4, Table 5; and Table 
4 and 5 in Paper I). 
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Table 5. Spectral fitting results - 111. Fits with complex models. 



Source No. Model^ 


r 




^edge/kT F— Statistic 


c.l.t 






(iO^i cm-2) 


(keV) 





96% 
>99.9% 
90% 
93% 
95% 
85% 
35% 
90% 



^ PL=SPL (or APL) + non-thermal model for the soft excess component; TM=SPL (or APL) + thermal model for the soft 
excess component; WA=SPL (or APL) + warm absorber features. * Confidence level with respect to model SPL (see Table 3) 
or model APL (see Table 4) using the F-statistic. 



43 WA 2.08I;;'* =Ng"' 0.55i;5'^° 4.7 

50 PL 1.66lE;'°^ 14-3l^i^ - 48.8 

58 PL 1.4llj;'^? 7.4l^-^ - 2.8 

63 PL 1.64lE;'g mtfi - 3.3 

72 PL 1.90f. 53.3l^g-^ - 4.6 

77 PL 1.33lE!'^^ 7.71^2^ - 2.1 

79 TM 2.06lJ;J^ 2.9l|? 0.43tg;|^ 0.4 

80 TM 1.90f. T-^til O-lS^o'io 2.6 



"dull" galaxies, i.e. source Nos. 3 and 41 (see Paper I), 
too. 

Before drawing conclusions from the results of the X- 
ray spectral analysis, we have checked out the possible 
presence of systematic trends due to the source position in 
the detector plane which could affect photon index and/or 
X-ray flux measurements. 

Fig. 3 shows the photon index obtained with the SPL 
model plotted against the hard X-ray flux: many sources 
have a r ~ 1.8-2 i.e. the typical value of unabsorbed AGNs 
(George et al. 2000; Nandra & Pounds 1994). This matches 
well with the optical identifications as well as with the re- 
sults obtained in the ASCA Large Sky Survey (Akiyama 
et al. 2000). On one hand, this plot shows that sources 
with flat slopes (F < 1.3) are present at various flux lev- 
els and no trend of F versus flux is evident in the data. 
On the other hand, despite the low number of objects at 
F2-10 > 10~i3 erg cm-2 g-^, very flat (i.e. F < 0.6) 
and inverted spectrum sources appear to be located in 



X 

0) 

TJ 2 


' "1 

: I 


1 1 ' 1 1 1 1 1 "1 


1 

i. - 


Photon 










'Oil 


1 1 





0.01 0,1 1 10 100 



Lum 2-10 10e(44) erg/s 

Fig. 4. Spectral index, computed assuming the best fit model 
for all optically identified sources in the sample (except for 
NGC 4291), as a function of hard X-ray absorption-corrected 
luminosity. 



3 FT 



I 




Fig. 3. Photon index obtained with the SPL model as a func- 
tion of the flux in the 2-10 keV band. 



the region of the lowest fluxes as predicted by the CXB 
synthesis model and recently observed in Chandra deep 
surveys (Tozzi et al. 2001; Stern et al. 2002). 

We have calculated the average SPL spectral indices in 
the 0.3-10 keV band for the BRIGHT and the FAINT sub- 
samples (see Sect. 2) and obtained (F) = 1.54 ± 0.03 and 
(F) = 1.56 ± 0.05, respectively. These results and their 
comparison with the values obtained from other hard X- 
ray surveys will be extensively discussed in Section 6.1. 

In Fig. 4 we plot the photon index as a function of 
the luminosity in the 2-10 keV band, as obtained apply- 
ing the best fit model for each optically identified source 
in the sample. There is no apparent trend for spectral 
variations as a function of luminosity. However, a large 
dispersion in the slope values is present: this could either 
be due to intrinsic differences (George et al. 2000) or to 
the contribution from additional spectral components (i.e. 
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Table 6. XMM-Newton properties of the 49 X-ray sources. 
Fluxes and luminosities are calculated for each source using 
the best-fit model listed in column 2 (see Table 5 for details). 



N 


Best-fit 
Model 


-'^O.S — 2 


-^2 — 10 




42 


SPL 


2.08 


3.79 




43 


WA 


1.21 


1.69 




44 


SPL 


2.73 


4.91 




45 


APL 


1.78 


4.38 




46 


SPL 


1.63 


2.91 




47 


SPL 


3.18 


5.01 




48 


APL 


3.07 


10.61 




49 


SPL 


1.73 


29.70 




50 


PL 


4.39 


28.1 




51 


APL 


3.83 


9.11 




52 


SPL 


1.63 


4.33 




53 


SPL 


4.46 


8.09 


22.7 


54 


SPL 


4.02 


4.87 




55 


SPL 


2.69 


3.31 




56 


APL 


0.52 


3.88 




57 


SPL 


1.71 


2.11 


10.8 


58 


PL 


0.95 


6.38 




59 


APL 


36.0 


76.0 


1.6 


60 


SPL 


10.64 


16.55 


2.2 


61 


SPL 


7.07 


4.08 


7.6 


62 


SPL 


8.97 


9.82 


5.1 


63 


PL 


1.18 


7.21 


5.4 


64 


SPL 


6.29 


4.13 


1.2 


65 


APL 


0.12 


3.98 


1.7 


66 


SPL 


3.63 


4.36 


27.1 


67 


SPL 


2.00 


3.40 


11.5 


68 


SPL 


5.92 


13.41 


3.6 


69 


APL 


0.04 


4.20 




70 


SPL 


2.49 


3.63 


5.9 


71 


SPL 


3.52 


5.57 


4.9 


72 


PL 


0.23 


6.76 




73 


SPL 


1.19 


1.14 




74 


APL 


0.07 


1.55 




75 


SPL 


1.73 


1.78 




76 


APL 


0.50 


1.18 




77 


PL 


0.74 


3.55 




78 


SPL 


1.51 


1.51 




79 


TM 


1.07 


1.03 




80 


TM 


0.53 


1.81 




81 


SPL 


1.43 


1.15 




82 


SPL 


2.64 


3.17 




83 


SPL 


1.17 


1.32 




84 


APL 


0.38 


2.16 




85 


SPL 


5.12 


2.49 


5.8 


86 


SPL 


3.04 


4.41 




87 


SPL 


3.27 


2.76 




88 


APL 


0.58 


1.73 




89 


APL 


1.92 


5.07 




90 


SPL 


1.96 


1.65 





in units of 10 erg cm ^ s ^ '''^ in units of 10** erg s 



reflection, soft-excess), unresolved here due to the Umited 
statistics. 



Redshift 

Fig. 5. Intrinsic 2-10 keV X-ray luminosity plotted against 
redshift for the identified sources in the sample (except for 
NGC 4291). The dashed line indicates the "completeness" flux 
limit of our survey (~ 5 x 10^^* erg cm~^ s^^) once all the 
FOVs are taken into account. 



Fig. 5 shows the luminosity-redshift relation for the 
optically identified sources in the sample: as expected, 
the minimum and the maximum hard X-ray luminosity 
in each redshift bin increases with the redshift owing to 
a selection effect (most luminous sources are detected far- 
thest). The present dataset is too small to disentangle any 
true luminosity dependence (as claimed in Barger et al. 
2002) from any mere selection effects. 

4.1. Multiwavelength properties of the sources 

A classical approach (Maccacaro et al. 1988) extensively 
used in the X-ray surveys to infer some information about 
the nature of unidentified sources as well as to test the 
reliability of the optical identifications themselves is the 
so-called "F_x/^bpt" diagnostic diagram. Since R mag- 
nitudes are available for the majority of the sources in 
the sample (Table 2), we adopt this band to build such 
a diagnostic diagram. The relation between optical mag- 
nitude and hard X-ray flux is plotted in Fig. 6. The flux 
in the R band is related to the R magnitude by the fol- 
lowing formula: logFji = —5.5 — 0.4 x R (Hornschemeier 
et al. 2001); hence, the relationship between optical flux 
and X-ray flux, obtained by Maccacaro et al. (1988) for 
the V magnitude, becomes the following: log{Fx / Fp) = 
logFx + 5.5 + R/2.5. 

The three dashed lines in Fig. 6 represent the X-ray to 
optical flux ratio of log{Fx /Fji) = 1,-1 and —2 (from top 
to bottom, respectively). They mark the regions occupied 
by different classes of X-ray sources as indicated in the 
Figure. A value comprised between —1 < log{Fx / Fn) < 
1 is typical of the "standard" luminous AGNs, both of 
Type 1 and Type 2. As expected on the basis of the op- 
tical identifications, most of the identified sources (indi- 
cated with triangles in the Figure) are properly located in 
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Fig. 6. Plot of optical 7? magnitude versus hard X-ray flux 
for those sources in the sample with available R band pho- 
tometry. The optically identifled (unidentifled) sources are in- 
dicated with triangles {pentagons) . The dashed diagonal lines 
have been calculated following Hornschemeier et al. (2001) and 
represent the constant flux ratios log{Fx / Fr) = 1 (upper), 
— 1 and —2 (lower). Labels show the typical {Fx/Fr) range 
observed for luminous AGN, starburst galaxies and normal 
galaxies in the local Universe. Peculiar sources (see text) are 
also indicated. 



the AGN locus. Interestingly, also most of the unidenti- 
fied sources (indicated with pentagons in the Figure) fall 
within the same region. This fact strengthens the results of 
our X-ray spectral analysis which indicated typical AGN- 
like spectra for all of them. These findings fully agree with 
results from hard X-ray surveys which found AGNs as the 
dominant population (e.g. Akiyama et al. 2000; Barger et 
al. 2002; Hasinger 2003). 

However, there are a few sources that have values of 
log{Fx / Ffj) outside the typical AGN locus. In particu- 
lar, there are 4 objects (i.e. Nos. 25, 38, 63 and 65) that 
are notably faint in the optical with respect to their X- 
ray flux values. They are of particular interest because re- 
sults from Chandra and XMM-Newton deep surveys (e.g. 
Alexander et al. 2001; Mainieri et al. 2002) suggest that 
this kind of sources are a mixture of "exotic" AGNs such 
as Type 2 QSOs, Extremely Red Objects (EROs, with 
R — K ^5), high redshift {z > 3) galaxies which host a 
dust-enshrouded AGN and, possibly, QSOs at z ^ 6. 

All the sources in the sample with log{Fx / Fj^) > 1, 
i.e. optically-weak, are found indeed significantly X-ray 
absorbed as shown in Tables 4 and 5, and, interestingly, 
they are also all RL QSOs: two of these are also classified 
as narrow line QSOs of the class of EROs (i.e. Nos. 63 and 
65), and the remaining two sources (i.e. Nos. 25 and 38) 
are optically classified as broad line QSOs. These latter 
show moderately flat X-ray spectra with best-fit photon 



indices F ~ 1.5-1.6 which are typically observed in core- 
dominated RL QSOs (Shastri et al. 1993; Sambruna et al. 
1999; Reeves & Turner 2000). Photometric information in 
the K band could be particularly useful to understand 
whether also these two QSOs could be classified as EROs 
as well. 

At the fluxes currently sampled by our survey, the 
population of objects with low ( ^ -1) and very low 
(< -2) values of log{Fx / Fr) are largely missed. Such 
kind of sources has been found to slowly emerge just be- 
low F2-10 ~ 10^^^ erg cm^^ s^^ and are thought to be 
mainly star-forming galaxies, LINERs and normal galax- 
ies (Barger et al. 2002). These classes of sources are com- 
monly detected in the soft X-ray band, where they are 
expected to be the dominat population at Fo.5-2 ~ 10~^^ 
erg cm~^ s~^ (Miyaji & Griffiths 2002). Three sources 
(Nos. 14, 73 and 85) in our sample show instead a ratio 
log{Fx / Ffi) < —1. Among them, the ffist two are opti- 
cally unidentified while the latter is a broad line QSO. 
The X-ray spectrum of this faint quasar (F2_io ~ 2.5 x 
10~^^ erg cm~^ s~^) is very steep (F ~ 2.6, e.g. Table 3), 
with the source counts fairly dominated by the soft X-ray 
photons. This is probably the reason why we measure a 
value of log{Fx / Fr) slightly lower than unity. 

As expected, the normal galaxy NGC 4291 (i.e. No. 
37) is the only object with a log{Fx / Fr) < —2. Besides 
normal galaxies, also nearby Compton-thick AGNs usu- 
ally exhibit such a low log{Fx / Fr) value (Comastri et al. 
2003). Accordingly, it is worth noting that the peculiar 
sources No. 3 (the "P3" galaxy, e.g. Fiore et al. 2000) and 
No. 41 (that are optically identified as "normal" galaxies, 
see Paper I) are found to lie in the region of the diagram 
typical of AGNs. This finding further supports the ob- 
scured AGN nature of these two objects. 

5. Hard X-ray selected QSOs 

A meaningful product of our wide-angle survey consists 
in the opportunity of investigating at faint flux levels the 
spectral properties of hard X-ray selected quasars (here- 
after HXSQs) and of comparing them with those derived 
from other studies of soft X-ray and/or optically selected 
QSOs. 

In order to extract QSOs from the list of sources with 
optical identification, we select those showing an intrinsic 
2-10 keV luminosity larger than 10^** erg/s. In doing so, 
we create a sample of 30 HXSQs that is one of the largest 



Table 7. Properties of hard X-ray selected QSOs. 



Objects 


N. 


O^Best~fit) 


(Fspl) 




AU QSOs 


30 


1.87±0.04 


1.66±0.03 


8.97 


QSOs at 2 <1 


20 


1.83±0.05 


1.52±0.04 


5.96 


QSOs at z>l 


10 


1.96±0.05'' 


1.96±0.05'' 


14.98 



" Luminosity in units of 10 erg s ^. 

'' Excluding source No.85: (rspL> = (rflest-/ri)=1.88±0.05. 
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sample of this kind available to date. Their redshifts range 
from z = 0.205 to z = 1.956 with a mean (z) ~ 0.7; the vast 
majority of the objetcs are broad lines AGNs with the re- 
markable exceptions of source Nos. 63 and 65, which show 
an optical Type 2 classification (Table 2; Mainieri et al. 
2002). The fraction of these HXSQs that turn out to be 
radio loud (i.e. with aoR > 0.3) is 20% (i.e. 6 out of 30 
sources) but, keeping in mind the incompleteness affecting 
the optical and radio coverage of our sample (see Table 2), 
this fraction should be considered only as a lower limit. 
Furthermore, we are able to split the sample of HXSQs in 
two subsamples having 10 and 20 objects each and red- 
shift larger and lower than 1, respectively. The average 
properties drawn from the spectral analysis of these two 
samples of HXSQs (total and subsamples) are reported 
in Table 7, while the photon indices obtained with the 
SPL and with the best fit model found in each individual 
HXSQ are displayed in Fig. 7 and Fig. 8, respectively. 

When the SPL model is applied we find a (F) = 1.66 
±0.03 which rises to (F) = 1.87±0.04 using the best fit 
model instead. This effect can be easily explained by the 
presence of intrinsic absorption which suppresses the soft 
portion of the X-ray emission. In particular, only HXSQs 
at z < 1, either radio-quiet (RQ) or radio-loud (RL), show 
this excess of absorption. Because of the limited number 
of HXSQs with z > 1, the corresponding value of (F) ap- 
pears to be biased by the very steep index of LBQS 2212- 
1759 (source No. 85). This source yields indeed a F 
2.6, similar to what found by Brinkmann et al. (2003) in 
another bright RQ QSO: such a steep slope is remines- 
cent of Narrow line Seyfert 1 galaxies, which usually show 
a strong soft excess component. However, any such soft 
excess component would be redshifted almost out of the 
EPIC energy range given the source redshift of z = 1.159. 
In any case, a more detailed modeling of such soft excess 
component is not allowed due to the low photon statistics. 
Excluding source No. 85 from the sample, we obtain in- 
deed an average index of (F) = 1.88±0.05, consistent with 
the average index found for the HXSQs at z < 1. We can 
therefore conclude that the average X-ray spectral slope 
of HXSQs resulting from the present work is almost the 
same from z « to z w 2 indicating no spectral evolu- 
tion, i.e. no obvious variation of F along the redshift (see 
Fig. 8). 

5.1. Comparison with previous works 

The value of F w 1.8-1.9 is similar to the mean slope of 
unobscured local Seyfert-like AGNs (i.e. F = 1.86±0.05, 
Nandra et al. 1997; Perola et al. 2002; Mahzia et al. 2003) 
and it matches well with previous X-ray studies of quasars 
carried out with different telescopes (Comastri et al. 1992; 
Lawson & Turner 1997). In particular. Reeves & Turner 
(2000) analysed a sample of 62 QSOs with ASCA, cov- 
ering a redshift range from 0.06 to 4.3, and they also re- 
ported a comparable mean photon index (F) — 1.76±0.04. 
They also claimed a difference between the slopes of RQ 
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Fig. 7. Plot of the full-band (0.3-10 keV) spectral index vs. 
redshift using SPL model for all the hard X-ray selected QSOs 
(HXSQs). The filled squares and the void stars are the "radio- 
quiet" HXSQs and the "radio-loud" HXSQs, respectively. It 
is worth stressing the incomplete radio coverage of our X-ray 
sample, so we include in "radio-quiet" HXSQs also some ob- 
jects with unknown radio properties. 
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Fig. 8. Same as in Figure 7, but with spectral indices derived 
using the best fit model for each source. 



QSOs (F - 1.9) and RL QSOs (F ~ 1.6), with the lat- 
ter showing a flatter average photon index. This issue 
has been extensively discussed in Sambruna, Eracleous & 
Mushotzky (1999, hereafter SEM99): these authors con- 
cluded that since lobe-dominated RL QSOs have the same 
intrinsic photon indices of RQ QSOs, the observed harder 
X-ray spectra of core-dominated RL QSOs being likely 
due to the presence of a beamed extra-continuum compo- 
nent originating in the inner parts of the radio jet. For the 
six RL HXSQs in our sample we find a (F) = 1.66±0.17 
which is slightly harder but still consistent with the typical 
value measured for RQ QSOs. However, it is worth not- 
ing that source Nos. 25 and 38, two bright objects with 
high-counting statistics, show flat photon indices (F ~ 
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1.3 and 1.5, respectively) similar to those usually found in 
core-dominated RL QSOs; however in these two sources 
we also find that the flat spectrum is due to the presence 
of intrinsic absorption. 

Our data also confirm that RL QSOs are characterized 
by absorption in excess to the Galactic value (Elvis et al. 
1994): 4 out of our 6 RL objects are indeed obscured by 
column densities of « 10^^"^^ cm~^. It has been suggested 
in other works (i.e. Cappi et al. 1997; SEM99), also on 
the basis of variability arguments, that such an absorber 
could be located in the inner nuclear regions of RL QSOs 
rather than associated with either the host galaxy or a sur- 
rounding cluster. Moreover, similar to our flnding, SEM99 
found that ^ 50% of broad line RL QSOs in their sam- 
ple suffered from significant intrinsic neutral absorption. 
These authors considered this fact at odds with simple 
orientation-based AGN unification models and suggested 
the possibility of different physical conditions of the gas 
around the central engine in RL QSOs, being colder in 
these objects than in RQ ones (George et al. 2000). 

5.2. Spectral evolution 

Concerning the evolution of the mean spectral shape, the 
value assessed for HXSQs by our analysis provides further 
support to the recent findings of Vignali et al. (2003). 
In fact, these authors invoke a universal X-ray emission 
mechanism for QSOs, i.e. indipendent from cosmic time 
and luminosity, on the basis of the analysis of nine high 
redshift (z > 4) QSOs observed with Chandra. By a 
"stacked" spectral fitting, they found an average spectral 
slope of 1.98 ± 0.16 in the 2-30 keV rest-frame band (but 
see also Bechtold et al. 2002). 

We report in Figure 9 a large compilation from litera- 
ture of photon indices obtained from the spectral analysis 
of radio-quiet QSOs up to z « 5.2 observed with different 
X-ray telescopes that includes also the recent Chandra re- 
sults mentioned above together with those found in the 
present work. No evident trend of the spectral slope ap- 
pears to emerge along with the redshift from this plot, 
confirming that the accretion mechanism in RQ QSOs is 
the same at any redshift sampled to date. 



6. Observational constraints for CXB synthesis 
models 

Two main observational constraints on the predictions of 
synthesis models can be extracted from the present work: 
(1) the average slope of the sources making the CXB at 
intermediate flux levels and (2) the corresponding ratio of 
absorbed-to-unabsorbed sources. Both these parameters 
are essential in the definition of the average spectrum at 
a given flux level in order to estimate the relative contri- 
butions of individual classes of objects to the CXB. 
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Fig. 9. Photon index versus redshift for a collection of RQ 
QSOs observed with different X-ray telescopes. Data points are 
taken, respectively, from: ASCA {asterisks, Reeves & Turner 
2000 and George et al. 2000), {void exagons, Vignali et al. 
1999); BeppoSAX {void triangles, Mineo et al. 2000). Data 
from the present work are indicated with filled exagons: we 
have excluded the QSOs from the PKS 0312-770 field be- 
cause of this field has no NVSS radio coverage. The void 
pentagons indicate XMMU J105144.6+572807 and XMMU 
J105125. 3-1-573048 (both at z ~ 3.4) which are the only RQ 
QSOs at z > 3 from the XMM-Newton long-exposure of the 
Lockman Hole (Mainieri et al. 2002). The void circles repre- 
sent two QSOs at z = 2.64 and z = 4.1 observed with XMM- 
Newton (Ferrero & Brinkmann 2003). The other three points 
at the highest redshifts correspond, respectively, to: APM 
08279-1-5222 {void square; z = 3.91, e.g. Hasinger et al. 2002), 
a joint spectrum of 9 Chandra QSOs with 4.09 ^ z ^ 4.51 
plotted here at their average redshift {filled circle, Vignali et 
al. 2003) and CXO J123647.9+620941 {void star ; z = 5.186, 
e.g. Vignah et al. 2002). 



6.1. Average slope at faint fluxes 

The average photon index calculated over the 0.3-10 keV 
band with the SPL model is (F) = 1.59 ± 0.02. As ex- 
pected, once the best fit model of each source is assumed^, 
the resulting value of the average slope becomes steeper 
i.e. with a (F) = 1.80 ± 0.04. This value agrees with 
the typical one of unabsorbed AGNs found in previous 
works (e.g. Lawson & Turner 1997; Reeves & Turner 2000; 
Malizia et al. 2003) and it is in agreement with the clas- 
sification as broad line AGNs for most of the identified 
sources (see Table 2). 

By dividing the whole sample into the BRIGHT and 
the FAINT subsamples, both of which are complete at 
their flux limit {F2-W 5 x lO^^'' erg cm^^ s^^ and 

® We exclude from this calculation NGC4291 (No. 37) as its 
spectrum is better described by a thermal model, in agreement 
with its optical classification as a normal galaxy (see Paper I). 
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Turner et al. (1997), i.e. a partial covering model with F — 
1.9, Nh — 10^^ cm^^ and a hard X-ray (deabsorbed) lu- 
minosity of 10'*'^ erg s~^. We have further assumed z = 0.7 
for the redshift value. The flux value calculated with these 
assumptions results to be equal to F2-10 « 4 x 10^^^ erg 
cm~^ s~^, i.e. well in the range where the spectral flat- 
tening has indeed been observed (Fig. 10). In addition, it 
is worth stressing that our simulations'' have shown that 
such a source would show a flat continuum if it were fitted 
with SPL model, i.e. as observed in an heavily absorbed 
object. 

If the bulk of the hard CXB originates from a large 
population of Compton-thin low-redshift Seyfert 2-like 
galaxies, they could indeed account for the flattening of 
the average photon index observed by Chandra towards 
very faint flux levels. 



Fig. 10. Average spectral index from a power model model 
(SPL) as a function of the hard X-ray flux. Data are taken 
from: Delia Ceca et al. (1999; exagon), Ueda et al. (1999; pen- 
tagon), Tozzi et al. (2001; stars). The triangle (square) repre- 
sents the value obtained in the present work for the BRIGHT 
(FAINT) subsample. 



-F2-10 XL 2 X 10 erg cm ^ s ^) we calculate in each 
case the average slope with the SPL model which turns out 
to be (F) = 1.54 ± 0.03 and (F) = 1.56 ± 0.05, respectively. 

In Fig. 10 the average photon indices of the BRIGHT 
and FAINT samples are compared with those obtained 
by the most popular hard X-ray surveys carried out so 
far. It is worth noting that our data resulted from the 
spectral analysis of individual sources, whereas all other 
values were derived by the stacking technique. A good 
agreement between our values and those obtained in other 
works is evident from this plot. 

At F2-10 ^ 1 X lO"^'' erg cm~^ s~^ the average slope 
appears still steeper than the integrated spectrum of the 
CXB, which has F = 1.4, thus suggesting that the bulk 
of the flat spectrum (i.e. absorbed) sources has still to 
emerge at these flux levels. The measurements obtained by 
the Chandra Deep Field South survey (Tozzi et al. 2001; 
indicated as stars in Fig. 10) confirm indeed this sugges- 
tion, revealing a progressive and significative fiattening of 
the average spectral index below F2_io ^ 10 x 10^^** erg 
cm~^ s~^, which is able to solve the "spectral paradox". 

As mentioned in Sect. 1, optical follow-up observa- 
tions of X-ray deep surveys have recently pointed out that 
the bulk of the hard CXB is produced by a large number 
of narrow line active galaxies at z 0.7 having Seyfert- 
like luminosities. Consequently, we have performed a test 
to show how the spectral hardening occuring at i^2-io 
10~^^ erg cm^^ s^^ could be qualitatively explained by 
these sources. Accordingly, we have assumed a typical X- 
ray spectral template of a Seyfert 2 galaxy as suggested in 



6.2. Fraction of absorbed sources 

Fig. 11 shows the fraction of absorbed (i.e. with an 
Nh ^ 10^^ cm^^ in excess to the Galactic value) to un- 
absorbed sources that we find in the BRIGHT and FAINT 
subsamples together with those obtained by other hard X- 
ray surveys with different flux limits in the flux range from 
F2-10 ~ 10~^^ down to -F2-10 ~ 10"^"* erg cm~^ s~^. 

Interestingly, it appears from this Figure that the frac- 
tion of absorbed sources remains almost the same 30%) 
in this wide range of hard X-ray fluxes. 

It is worth noting that in order to calculate the frac- 
tions for both samples reported in this Figure we assume 
F = 1.9 for all the sources and z — 1 for the unidentified 
ones to overcome some biases which may affect our esti- 
mates, as mentioned in Sect. 3.1.1. The values reported in 
Fig. 11 can be therefore considered as conservative esti- 
mates of the fraction of absorbed sources in both samples. 

6.2.1. Comparison with theoretical predictions 

This finding is fairly unexpected since all synthesis mod- 
els of the CXB (e.g. Gflfi et al. 2001; Wilman & Fabian 
1999; Comastri et al. 2001, hereafter COl) predict that 
as fainter fiuxes are considered, more absorbed sources 
should be found. 

In particular, in Figure 12 the theoretical predictions of 
the CXB synthesis model of COl are plotted together with 
the fraction of absorbed sources in the 2-10 keV band 
found in our spectral survey. For a better comparison we 
have splitted our measurement regarding the BRIGHT 
sample into three points. Moreover we have added the 
point relative to 9 x lO'^^ < F2-10 ^ 2 x 10"" erg 
cm^^ s"^ calculated taking into account all the sources 
selected in the field of LBQS 2212-1759 (this field has in- 
deed the longest exposure time amongst the observations 
considered here, i.e. 80.5 ks, see Table 1) and those in 

^ Assuming an input Seyfert 2-Iike spectrum as reported 
above and a 300 ks PN exposure, we found a photon index 
r = 0.92lo;43 in the case of a fit with the SPL model. 
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Fig. 11. Fraction of absorbed sources (i.e. witli Nh ^ 10^^ 
cm~^) taken from different iiard X-ray surveys. The filled 
square and the filled circle indicate the values found from our 
BRIGHT (60% optically identified, hereafter ID) and FAINT 
(~ 35% ID) samples, respectively. Other symbols represent the 
values from: Piccinotti sample (Gomastri 2000, but also in- 
cluding 4 BL Lacs and M82 present in the original sample of 
Piccinotti et al. 1982; filled pentagon; 100% ID), BeppoSAX 2- 
10 keV survey (Giommi et al. 2000); void star), ASCA Large 
Area Survey (Akiyama et al. 2000; void triangle; 97% ID), and 
XMM-Newton survey in the Lockman Hole (Mainieri et al. 
2002, void exagon; ~ 90% ID). The horizontal error bars as- 
sociated to the points represent the flux interval between the 
completeness limit and the highest flux sampled in each survey; 
the vertical error bars indicate instead the Poissonian errors for 
each fractional value. 



the 100 ks exposure of the Lockman Hole presented in 
Mainieri et aL (2002). 

This figure shows that the fraction of absorbed sources 
predicted by the theory is a factor of ~ 2 larger (with a « 
4.5cr significance) than our observational data at compa- 
rable fluxes. We thus confirm and extend even to fainter 
fluxes the data-to-model mismatch found in Paper I. Our 
data suggest that unabsorbed (i.e. with Nh ^ 10^^ cm~^) 
objects still largely dominate the source counts at -F2-10 
^ 10~^^ erg cm~^ s~^, at odds with the theoretical ex- 
pectations of about 50-60% of absorbed sources. 

Despite this mismatch, the overall trend predicted by 
the model follows the observational data as if only the 
normalization was wrong. 

Although very preliminary, also the results of the in- 
dividual spectral analysis of the brightest X-ray sources 
in the HDFN (Bauer et al. 2003) seem to suggest a 
Nh distribution peaked towards low column density val- 
ues with only few sources (< 9%) obscured by Nh ^ 10^^ 
cm~^ whereas the synthesis models predict a value of 
~ 20%. Other works based on hardness ratios analy- 
ses further confirm such a scarsity of obscured objects 
above F2-10 ~ lO^^** erg cm~^ (e.g. Baldi et al. 
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Fig. 12. Fraction of absorbed sources (i.e. with Nh 



10^ 



cm~ ) from our analysis vs. 2-10 keV flux compared to theo- 
retical prediction of Gomastri et al. (2001). For the calculation 
of our points we assumed a conservative spectral slope of r= 1 .9 
for all sources and z = 1 whenever the X-ray source was not 
optically identified. The value for the point at faintest flux has 
been calculated using sources selected in the deepest exposure 
of ours (i.e. the LBQS 2212-1759 filed, 80.5 ks) together with 
sources in the 100 ks exposure of the Lockman Hole analysed 
by Mainieri et al. (2002). 



2002; Akiyama, Ueda & Otha 2002): in particular the lat- 
ter authors reported a fraction of absorbed source with 
Nh ^ 10^2 cm-2 and L2-10 > 10'*'* erg s"! of - 15%, 
i.e. a factor of 3 lower than expected. 

Furthermore, the average photon index values at 
-F2-10 ~ 10"-^^ erg s~^ reported in the present as well 
as in other works (Sect. 6.1) result to be steeper than the 
slope of the integrated CXB spectrum: this result is indica- 
tive of the fact that absorbed objects are not yet present 
in large quantities at these fluxes, as suggested instead by 
theoretical models. 

Possible explanations for this observational mismatch 
will be extensively discussed in Sect. 7. 



6.3. Comparison with the results from the IMs HDFN 
survey 

Since this data/model mismatch on the fraction of ab- 
sorbed sources has been claimed in Paper I and in the 
present paper for the first time, it requires further inves- 
tigations before validation. We have thus used the data 
from the Chandra IMs Hubble Deep Field North survey 
published in Brandt et al. (2001) to estimate the fraction 
of absorbed versus unabsorbed objects at our flux levels 
and lower. 

To this aim, from the entire X-ray source catalog only 
the sources detected in the 2-8 keV band (265) have been 
selected. Their fluxes span a very large range i.e. 1.4 x 
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10-16 < F2-8 ;S 1.2 X 10"" ergcni-2 s-i. Ofthese265 
sources, Barger et al. (2002) presented the optical identi- 
fication for 118 (i.e. ~ 45%). 

Since Brandt et al. (2001) used the source fluxes cor- 
rected for vignetting, we calculate source-by-source the 
"flux ratio" FR rather than the commonly-used "hard- 
ness ratio". This ratio is defined as follows: 

FR = {F2-S - ^0.5-2)/(-F2-8 + Fo.5-2) 

Similarly to the "hardness ratio" , this quantity is indica- 
tive of the "flatness" of an X-ray spectrum. Accordingly, 
sources with FR = 1 represent those without a positive 
detection in the soft X-ray band but seen in the hard X- 
ray band. We then compare the FR value of every source 
with that (hereafter indicated with FR^) obtained^ for 
an object with F = 1.7, Nh = lO'^^ cm~^ placed at the 
redshift of that source (or z = 1 for an optically unidenti- 
fied one). Accordingly, a source that shows FR > FRz is 
X-ray obscured with an Nh > 10^^ cm~^. 

Following this procedure, it is possible to estimate the 
fraction of absorbed sources in the Chandra HDFN IMs 
exposure at different hard X-ray flux levels. These results 
are displayed in Figure 13 together with those derived 
from our analysis. 

It appears from Figure 13 that at hard X-ray fluxes 
from « 10""'^'* erg cm~^ s"^ down to « 10"^^ erg cm~^ 
s^i the fraction calculated in the HDFN data matches 
very well with the model predictions. On the other hand, 
the two points at F2-S > 10"^^ erg cm^^ s^^ are consis- 
tent both with the theoretical values and our estimates. 
However, it is worth noting the large errors which affect 
the Chandra measurements: in fact, due to their lower 
number densities per deg^, the number of objects detected 
at the brightest fluxes by this 1 Ms exposure is relatively 
small, with only 34 X-ray sources having F2-S > 10"^^ 
erg cm~^ s^^. Furthermore, some flat spectra could be 
due to warm absorber features rather than cold absorp- 
tion material and unfortunately hardness ratio analysis 
does not allow to discriminate between the two. 

The analysis of this ultra-deep Chandra obser- 
vation therefore confirms only partially our finding. 
Unfortunately, it does not provide an efficient tool to dis- 
entangle the mismatch reported in Section 6.2. 

This result is however very useful as it provides the 
most accurate estimate never published elsewhere for the 
fraction of absorbed sources at flux levels much fainter 
than those sampled in our survey ( ^ 10^^^ erg cm^^ s^^). 

7. How to explain this data/model mismatch? 

7.1. Observational biases 

There are three main observational biases which may af- 
fect our estimate of the fraction of obscured sources, i.e.: 



® These values of FRz are determined using the A02 version 
of PIMMS (Mukai 2000). For example, for a source with V = 



IIDI''N IMS (i'-l 7, z-rcai or z-l] 
model prodicLion (CoriiasLi'i 01) 



|(8-15)e-14 

h), I : _ 
I (a 5)e 14 
^ (0.9-2)e- 14 
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Flux (2-8 keV) cgs 

Fig. 13. Ratios of absorbed (Nh W'^^ cm~^) to unabsorbed 
hard X-ray selected sources in our sample {filled symbols) and 
in the HDFN IMs exposures {void pentagons) as a function 
of the 2-8 keV flux compared to model predictions {dashed 
line; Comastri et al. 2001). Points from our analysis are the 
same as those plotted in Fig. 12, but at flux values smaller 
than a factor of 1.2 (F2-10/F2-8 ~ 1.2 with F = 1.7) in order 
to provide a better comparison with HDFN data. The val- 
ues from the HDFN has been derived taking into account X- 
ray and optical data from Brandt et al. (2001) and Barger et 
al. (2002), respectively. We obtained the fractions of absorbed 
sources in each flux bin by means the "flux-ratio", i.e.: FR 
= (-f2-8—-Fo. 5-2)7(^2-8-1-^0. 5-2), analysis. We consider "ab- 
sorbed" a source with an FR value larger than that calculated 
assuming F — 1.7, Nh = 10^^ cm~^ and the redshift proper 
of the source or 2: = 1 if the X-ray source is not optically 
identified (see text for further details). 



1.7, Nh 
is 0.41 



10 cm and z =1, the corresponding FRz value 



Poor statistics from the faintest sources. This fact may 
have lead in some cases to a possible underestimate of 
the real Nh value because of the possible presence of 
an unresolved soft-excess or reflection component. 
However, we consider that this bias is in large part 
taken into account by fixing the photon index to F = 
1.9 before re-estimating the amount of cold absorption 
column density (see Sect. 3.1.1). Moreover, as reported 
in Sect. 4.2.2 of Paper I, we performed simulations in 
order to rule out the possibility that we miss some 
obscured sources at fluxes above the completeness limit 
of our survey. 

Incompleteness of the optical identifications. The col- 
umn density measured in the optically unindentified 
sources would be certainly underestimated if they were 
placed at z = (see Sect. 3.1.1). Nevertheless, we es- 
timate that we overcome this bias by assuming z — 1 
for such sources, i.e. by placing them just above the 
peak of the newly observed redshift distribution (e.g. 
Hasinger 2003; Barger et al. 2002) in order to provide 
a more conservative estimate of their Nh values. 
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3. Selection effect due to the PN effective area. The pre- 
dictions shown in Fig. 12 and Fig. 13 have been cal- 
culated assuming a flat-response detector (COl). This 
is not the case of XMM-Newton, whose effective area, 
although the best to date, is affected by a degrada- 
tion by about a factor of two between 4 and 9 kcV. 
Hence, a convolution of the models with the response 
of EPIC would be very important in order to allow a 
proper comparison between the observational data and 
the theoretical expectations^. 

Interestingly, if we rescrict the analysis to the 5-10 

keV energy band and compare the fraction of absorbed 
objects resulting from the complete subsample of 22 
sources with F^_io ^ 5 x 10~^^ erg cm~^ s~^ with 
the model prediction in this band, we still notice that 
our value is ~ 30%, i.e. a factor of 2 lower than ex- 
pected by COl (~ 60%). 

7.2. Revision of some model assumptions 

We address here the possibility that one (or more) of the 
assumptions usually included in synthesis models of the 
CXB need to be revised. The results presented here and 
the mismatch emerging in the optical follow-up of the 
deep X-ray surveys (Hasinger 2003) between the predicted 
and the observed redshift distribution of X-ray sources 
indeed suggest that some inputs of the CXB theoretical 
models should be updated. 

The most critical assumptions usually made in theo- 
retical works are three: the spectral template adopted for 
the X-ray sources, the distribution of absorption column 
densities and the slope, form and evolution of the X-ray 
luminosity function (XLF), which is almost unknown for 
Type 2 objects. 

1. The X-ray spectral properties. Although the overall X 
ray spectral shape of AGN is thought to be known 
(Mushotzky et al. 1993, Nandra & Pounds 1994), some 
details arc still uncertain. 

For example, the average value of the high-energy cut- 
off (due to thermal Comptonization) is poorly con- 
strained. Most theoretical models assume Ecutoff = 
320 keV (COl, GiUi et al. 1999), but this is certainly a 
rough approximation. On the basis of BeppoSAX ob- 
servations Matt (2000), Perola et al. (2002) and, re- 
cently, Malizia et al. (2003) clearly showed that the 
average cutoff value is more around 120-180 keV. 
It is worth noting that the high-energy cutoff is one 
of the key parameters in reproducing the peculiar 
"bump-like" shape of the CXB. The ubiquitous pres- 
ence of an upward curvature (the so-called soft-excess) 

^ Interestingly, Gilli et al. (2001) showed that the fraction 
of obscured AGNs (i.e. with Nh > 10^^ cm^^) predicted by 
their model at ^2-10 ~ lO^^"* erg cm^^ s^^ decreases by 20% 
when folding the model predictions with the ASCA SIS effec- 
tive area. If the fraction of obscured AGNs predicted by COl 
is decreased by 20% at ^2-10 w 10~^^ erg cm~^ the dis- 
crepancy observed in Fig. 12 would be consequently reduced. 



rising steeply below 2 keV in the spectrum of unab- 
sorbed AGNs is also uncertain (Matt 2000), especially 
in the case of luminous objects (Reeves & Turner 2000, 
George et al. 2000). 

All the above spectral properties are still to be well 

constrained for the brightest sources in the local 
Universe and, mostly, there is a evident lack of infor- 
mation about their possible evolution along z and/or 
luminosity. 

2. The Nu distribution. GiUi et al. (1999) and Gilh et 
al. (2001) assumed in their models a "fixed" distri- 
bution of the absorbing column densities for Type 2 
objects. They used the distribution found by Risaliti 
et al. (1999) for a sample of nearby optically-selected 
Seyfert 2 galaxies. There are three main causes of un- 
certainty in their assumption: (z) the observed distri- 
bution is based on a limited luminosity range; (m) it 
takes into account only "optical" Type 2 AGNs i.e. ob- 
jects with narrow emission lines in their optical spectra 
(but, as extensively discussed in Paper I, a mismatch 
between optical and X-ray classification of AGN has 
been widely observed) and, most important, (Hi) such 
distribution has been derived from sources in the local 
Universe (which are clearly not responsible for the bulk 
of the CXB). Consequently, the assumptions made in 
their model regarding the extension of the Nh distribu- 
tion to higher redshifts could be erroneous, thus intro- 
ducing an error in the final output of their estimation. 

3. The XLF of Type 2 AGNs. Recent results from opti- 
cal follow-ups of Chandra and XMM-Newton deep 
surveys yield a source redshift distribution that, if 
confirmed, would require a substantial revision of the 
commonly assumed XLFs (Hasinger 2003; Cowie et al. 
2003). 

The main uncertainty in all synthesis models of the 
CXB is indeed the XLF of obscured sources and its 
cosmological evolution, both of which are completely 
unknown. So far theoretical models have assumed the 
same XLF for Type 1 and Type 2 AGNs. However, on 
the basis of the above results (i.e. the bulk of the CXB 
originates from narrow line AGNs at z ^ 1), it has 
been suggested that the evolution properties of Type 
2 sources is likely to be different from that of Type Is. 
Moreover, the c;ommonly-used XLFs of Type 1 AGNs 
have been mainly derived from soft X-ray surveys 
and, in addition, their evolution with z is still debated 
(Comastri 2000 for a complete review). 
An attempt to adapt the theory to the new ob- 
servational evidences, has been recently done by 
Franceschini, Braito fc Fadda (2002, hereafter FBF02). 
These authors have proposed a model for the CXB 
where the obscured AGNs closely follow the evolution 
of strongly evolving infrared starburst galaxies (as re- 
cently found by ISO surveys), which evolve steeply up 
to ^; ~ 0.8 similarly to the X-ray sources detected in 
the deep surveys (Cowie et al. 2003). This approach is 
based on the fact that ~ 65% of the sources (mostly 
AGNs) detected in the Lockman Hole survey in the 5- 
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10 keV band have an IR counterpart at 15 fjtm.. FBF02 
therefore suggested two different evolutionary patterns 
for Type 1 and Type 2 AGNs: the former evolve as 
given by optical and soft X-ray surveys, whereas the 
latter evolve faster, as found in mid-IR surveys. 
However, this scenario seems to be more complicated 
than described by FBF02 (and subsequently refined 
by Gandhi & Fabian 2003): in fact, GilU (2003) has 
shown that such a "starburst-like" evolution for Type 
2 objects largely overestimates the ratio of obscured- 
to-unobscured objects at ^; ^ 1. 
4. Space density evolution of different AGN types. 
Another hypothesis is that classes of AGN with dif- 
ferent luminosity may evolve differently, contrary to 
what generally assumed in theoretical models, where 
the evolution is modeled only taking into account un- 
absorbed QSOs. 

In particular, Cowic ct al. (2003) and Bargcr et 
al. (2003) have recently suggested on the basis of 
Chandra data that the space density evolution for 
Seyfert-like objects could be nearly constant (or, at 
most, slightly dechning above z > 0.7) up to z ~ 2.5, 
while that of high-luminosity QSO-like AGNs shall in- 
crease rapidly from z ^ up to 2 3 in agreement 
with the evolution of optically-selected samples. 

The need for more accurate XLFs both for obscured and 
unobscured sources clearly emerges from this discussion. 
This information is indeed a crucial input for all models 
aimed at achieving a correct and complete interpretation 
of the CXB phenomenon. 

To date, it appears difficult to predict how all this new 
information obtained by Chandra and XMM-Newton will 
affect synthesis models of the CXB. It is beyond the goal 
of this paper to estimate what type of changes these new 
observational constraints (including our new measurement 
of the fraction of absorbed versus unabsorbed objects) will 
put on these models, but it will be clearly an important 
development for future theoretical works. 

8. Conclusions 

This work provides the first step in the detailed study 
of the X- ray spectral properties of hard X-ray selected 
sources detected at faint fluxes (i.e. F2_io ^ 10"^'^ erg 
cm~^ s~^) near the knee of the LogiV Logs' distribution. 
These are the sources that most contribute to the CXB. 

Previously published works on this topic were based 
mainly on hardness ratio and/or stacked spectral anal- 
yses. Complementary to ultradeep pencil-beam sm-veys, 
our shallower survey addresses for the first time the anal- 
ysis of each individual spectrum. 

Results have been reported for a sample of 90 hard 
X-ray selected sources detected serendipitously in twelve 
EPIC fields. This is the largest ever made sample of this 
type. A detailed spectral analysis has been performed in 
order to measure source-by-source the 0.3-10 kcV contin- 
uum shape, the amount of cold (and, possibly, ionized) 



absorbing matter and the strength of other spectral fea- 
tures such as soft excess and warm absorber components. 

The most important results can be briefly outlined as 
follows: 

1. Fluxes in the 2-10 (0.5-2) keV band span from 1 (0.04) 

to 80 (70) X 10"" erg cm^^ g-i. About 40% of the X- 
ray sources are optically classifled from the literature. 
Most of them arc broad line AGNs with redshift in 
the range 0.1 ^ z ^ 2. The high luminosities foimd 
(lO'*^ < L2-10 < 10^^ erg s-i) match well with these 
identifications except for two "optically dull" galaxies 
(see also Paper I). 

2. Using a simple power law model with Galactic ab- 
sorption we obtain (F) = 1.59±0.02. Considering only 
sources in the BRIGHT and FAINT subsamples we 
find (F) = 1.53±0.03 and (F) = 1.56±0.05, respec- 
tively. Both these values are in fairly good agreement 
with stacked spectral analysis obtained from ASCA 
and Chandra hard X-ray surveys. 

3. 65% of the sources are well fitted with the SPL model; 
their average spectrum provides a photon index F ~ 
1.7-^2.0, i.e. typical of unabsorbed AGNs. 

4. 30% of the sources require a column density larger than 
the Galactic value with Nh ranging from « lO^'^ to 
PS 10^"^ cm~^ (Sect. 4). In particular, two narrow line 
AGNs turn out to be Type 2 QSOs since they are char- 
acterized by high Imninosities (L2-1Q > 10^^ erg s"""^) 
and high column densities (Nh > 10^^ cm^^). 

5. The mean slope of hard X-ray selected QSOs in our 
sample remains nearly constant ((F) « 1.8 1.9) be- 
tween z ~ and ^ 2 (Sect. 5). By combining this 
result with other recent works on high z QSOs, we 
strengthen the suggestion that QSOs do not exhibit 
any spectral evolution and, hence, the type of accre- 
tion in these objects should be similar up to z « 5. 

6. While from the analysis of the sources detected in the 
IMs HDFN survey at faint fluxes (^2-10 < 10" erg 
cm~^ s~^) the observed fraction of absorbed sources 
(Nh ~ 10^^ cm~^) is consistent with the theoretical 
predictions, at the brighter fluxes (i^2- 10 ^ 10~^^ erg 
cm~^ s~^) considered in our survey it appears to be a 
factor ~ 2 lower (with a « 4.5c7 significance) than pre- 
dicted by synthesis models of the CXB. This confirms 
and extends our previous results obtained in Paper I. 
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